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ABSTRACT 



Context. LS I +61° 303 is one of the most interesting high-mass X-ray binaries owing to its spatially resolved radio emission and its 
TeV emission, generally attributed to non-thermal particles in an accretion-powered relativistic jet or in the termination shock of the 
relativistic wind of a young pulsar. Also, the nature of the compact object is still debated. Only LS 5039 and PSRB1259-63 (which 
hosts a non-accreting millisecond pulsar) have similar characteristics. 

Aims. We study the X-ray emission from LS I +61° 303 covering both short-term and orbital variability. We also investigate the source 
spectral properties in the soft X-ray (0.3-10 keV) energy range. 

Methods. Twenty-five snapshot observations of LS I +61° 303 were collected in 2006 with the XRT instrument on-board the Swift 
satellite over a period of four months, corresponding to about five orbital cycles. Since individual data sets have too few counts for a 
meaningful spectral analysis, we extracted a cumulative spectrum. 

Results. The count rate folded at the orbital phase shows a clear modulation pattern at the 26.5 days period and suggests that the 
X-ray peak occurs around phase 0.65. Moreover, the X-ray emission appears to be variable on a timescale of ~1 ks. The cumulative 
spectrum is well described by an absorbed power-law model, with hydrogen column density N H = (5.7 ± 0.3) x 10 21 cirr 2 and photon 
index T = 1.78 + 0.05. No accretion disk signatures, such as an iron line, are found in the spectrum. 

Key words. X-rays: individual: LS I +61° 303- X-rays: binaries 



1. Introduction 

LSI +61° 303 is a peculiar binary system characterized by 
variable radio emission, long known to be modulated with a 
26.4960 ± 0.00 28 days orbital period dGregorv & Taylor! 1 19781; 
lGregoryil2002h . The eccentricity of the orbit is 0.72 ± 0.15, 
and the periastron passage occurs at phase <p = 0.23 + 0.02 
dCasares et alj|200"5l) . The phase of the radio maxim um shifts 
betwe en 0.45 < <f> < 0.9 with a 1667 + 8 days period (Gregory 
|2002|) . pointing to a second modulation. H 1 observations of the 
radio source give a distance of 2.0 + 0.2 kpc dFrail & Hiell ming 
1 1 99 lb . implying an X-ray luminosity of ~10 33 erg s _1 . Optical 
studies have established the primary star to be a Be surrounded 
by a small and inhomogeneou s disk whose optical and Ha 
emission is also modulated dH utchings & Crampton 1981; 
Men delson & Mazehl [19891 iParedes et al.lll994t IZamanov et alj 



as a gamma-ray emitter ( Alber t~et al.l 120061) . This was also 
confirmed in the hard X-ray energies (up to 100 keV) by 
INTEGRAL observation of variable e mission, peaking a t phase 
~0.55 of the orbital m odulation dChernvakova et all 2006; 
iHermsen & K uiper 2007|). 

Although in soft X-rays there were no doubts from the 
start about the source identification (Bignami e t al J [19811) . the 
X-ray observations, aimed at studying the orbital modulation, 
proceeded in parallel with the gamma-ray ones, using all 
the X-ray instrument s active ov er the past 30 years: Einstein 
dBignami et all 1 198 lb . ROSA T dGoidoni & Mereghettil 119951; 



Tavlo r~et~aT] 



(Paredes et al 
1200 U iLeahv 




A SCA dLeahvet al.1 119971). RossiXTE 



(Chernyakova et al. 



1999). 



iHarrison et alj 20001; Greiner & Raul 
iGrundstrom et alj 120071). IN T EGRAL 
120061 IHermsen & Kuipeil 120071) . 

Chandra 



LSI +61° 303 has been attracting a lot of attention (and 
observing time) owing to its putative link with the gamma-ray 
sourc e 2CG 135+01, discovere d by COS-B ( Hermse n et alj 
119771; ISwanenburg etafl 1 198 It iBignami & Hermsenl 1 1983b 
30 years ago and later confirmed by CGRO/EGRET, as 
3EGJ0241+6103 dHartman et all 119991) . with observations 
spanning the period 1991-2000. In s pite of relentless efforts 
dTavani et alj [l996; Knif fen et al.l ll997). no conclusive evidence 
for the orbital modulation was found in the gamma-ray data, 
thus hampering the source identification. Recently the orbital 
modulation was seen by the MAGIC Cherenkov air-showers 
telescope in ultra high energy gamma rays (above 400 GeV), 
clinching the case for the identification of LSI +61° 303 
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BeppoSAX, XMM-N ewton dSidoli et all 120061) . and 
(IParedes et al.ll2007l) . 

Today, the behavior of LSI +61° 303 continues to defy 
easy characterization. Its X-ray emission ap pears to be variabl e 
over both long and short time scales (e.g.. ISidoli et al.l l2006). 
hampering the detection of an unambiguous modulation pattern 
when only sparse observations are available. Here we report the 
monitoring of the source soft X-ray flux over 5 adjacent orbital 
cycles, using t he XRT telescope onboard the NASA/UK/ASI 
Swift mission dGehrels et al.ll2004l) . 



2. Observations and analysis 

Twenty-five Swift observations of LS I +61° 303 have been per- 
formed with the XR T instrument in photon counting (PC) mode 
dBurrows etal.l2 005). The XRT uses a CCD detector sensitive to 



2 

Table 1. Observation log. 
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Sequence 


Start time (UT) 


End time (UT) 


Phase" 


Net exposure b 


Count rate c 




(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 




(s) 


(counts s -1 ) 


00036093001 


2006-09-02 08:09:31 


2006-09-02 11:46:56 


0.591-0.597 


4665 


0.271 


± 0.008 


00036093002 


2006-09-05 21:39:40 


2006-09-05 23:22:56 


0.726-0.728 


749 


0.22 


±0.02 


00036093003 


2006-09-11 00:57:04 


2006-09-11 23:30:56 


0.920-0.955 


4657 


0.120 


± 0.005 


00036093004 


2006-09-13 09:40:21 


2006-09-13 22:36:57 


0.009-0.029 


3598 


0.187 


± 0.007 


00036093005 


2006-09-15 07:56:34 


2006-09-15 13:11:57 


0.082-0.090 


3783 


0.223 


± 0.008 


00036093006 


2006-09-17 06:33:09 


2006-09-17 09:54:58 


0.155-0.160 


3559 


0.091 


± 0.005 


00036093007 


2006-09-19 09:53:05 


2006-09-19 14:57:57 


0.236-0.244 


2053 


0.117 


± 0.008 


00036093008 


2006-11-21 01:56:37 


2006-11-21 08:22:56 


0.601-0.611 


3024 


0.226 


± 0.009 


00036093009 


2006-11-22 03:28:09 


2006-11-22 08:27:56 


0.641-0.649 


3965 


0.219 


± 0.007 


00036093010 


2006-11-23 01:59:58 


2006-11-23 07:01:56 


0.676-0.684 


3092 


0.30 


±0.01 


00036093011 


2006-11-24 03:47:00 


2006-11-24 08:46:56 


0.717-0.725 


2976 


0.30 


±0.01 


00036093012 


2006-11-29 06:13:28 


2006-11-29 08:02:56 


0.909-0.912 


1126 


0.13 


±0.01 


00036093013 


2006-11-30 06:18:01 


2006-11-30 08:06:56 


0.947-0.950 


1900 


0.22 


±0.01 


00036093014 d 


2006-12-02 22:45:00 


2006-12-02 22:45:57 


0.049 


58 


0.24 


±0.07 


00036093016 


2006-12-05 03:42:40 


2006-12-05 13:23:58 


0.132-0.147 


1284 


0.12 


+ 0.01 


00036093017 


2006-12-07 03:47:46 


2006-12-07 16:51:57 


0.208-0.228 


2455 


0.065 


± 0.005 


00036093018 


2006-12-09 12:24:07 


2006-12-09 16:59:57 


0.297-0.304 


2145 


0.20 


±0.01 


onrnfsncnni o 

UUUJUUyJU 1 y 


?nnfi c\ci-^ 1 -or 

ZUUU -1Z.-11 UU.J1 .uo 


ZuuU -1Z -1 1 Uil.H'O.JO 


U.JJH — U.JJ / 


?no7 


0.076 


± 0.006 


00036093020 


2006-12-13 01:01:26 


2006-12-13 09:07:57 


0.430-0.442 


1126 


0.12 


±0.01 


00036093021 


2006-12-14 01:09:47 


2006-12-14 23:40:57 


0.468-0.503 


707 


0.14 


±0.01 


00036093022 


2006-12-16 01:21:46 


2006-12-16 14:14:56 


0.543-0.564 


1633 


0.19 


±0.01 


00036093023 


2006-12-18 06:15:42 


2006-12-18 17:36:58 


0.627-0.644 


2168 


0.26 


±0.01 


00036093024 


2006-12-20 19:19:35 


2006-12-20 19:39:56 


0.722-0.723 


622 


0.26 


±0.02 


00036093025 


2006-12-22 00:12:58 


2006-12-22 05:07:58 


0.768-0.776 


1393 


0.13 


±0.01 


00036093026 


2006-12-24 10:27:07 


2006-12-24 15:21:56 


0.860-0.867 


2764 


0.25 


±0.01 



a Following Gregory! j2002h we adopt as time of the zero phase JD = 2, 443, 366.775, with the orbital period P orb = 26.4960 ± 0.0028 days. The 
mean error on the orbital phase i s ±0.002 . The phase of the super-orbital modulation of the orbital phase and amplitude of the radio outbursts 
(with a period of 1667 ± 8 davs: lGregorvll2002h spans from -0.367 (observation 00036093001) to -0.435 (observation 00036093026). 

b The exposure time is spread over several snapshots (single continuous pointings at the target) during each observation. 

c Hereafter all errors are at 1 cr confidence level, unless otherwise specified. 

d This observation is listed for completeness, but will not be used for the analysis owing to its short integration time. 



photons with energies between 0.2 and 10 keV. Table Q] reports 
the log of the Swift/XRT observations used for this work. 

The data were processed with standard procedures us- 
ing the FTOOLS0 task xrtpipeline (version 0.10.6 under 
the Heasoft package 6.2.0). We selected events with grades 
0-12 and limited the analysis between 0.3-10 keV, where the 
PC response matrices are well calibrated. We extracted the 
source events from a circular region with a radius of 20 pixels 
(1 pixel =! 2.37"), whereas to account for the background we ex- 
tracted events within an annular source-free region centered on 
LSI +61° 303 and with radii 55 and 75 pixels. Since the data 
show a maximum count rate <0.3 counts s , no pile-up correc- 
tion was necessary. 

Inspection of the XRT light curves relative to each obser- 
vation shows evidence of moderate variability (up to a factor 
of ~4) on a timescale of -1 ks (see Fig.[TJi, an eff ect already 
obser ved in X-rays with various instruments (e.g., Sid oli et alj 
120061) . However, the statistics are too poor to assess spectral vari- 
ation as a function of the source count rate. 

Individual data sets have too few counts for a meaningful 
spectral analysis. Therefore, we extracted a cumulative spec- 
trum. This corresponds to a total exposure of 52.6 ks and con- 
tains about 10,000 net counts in the 0.3-10 keV range. The 
data were rebinned with a minimum of 20 counts per en- 
ergy bin to allow x 1 fitting. The ancillary response file was 



1 See http : //heasarc . gsf c . nasa . gov/ docs/software/ . 
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Fig. 1. Swift/XRT light curve of the first 14 ks of observation 
00036093003 (see Table[TJfor details). The bin size is 250 s and 
less than 0.001 counts s _1 are expected from background events. 
The short-term X-ray variability of LS I +61° 303 is apparent. 



generated with xrtmkarf, and it accounts for different ex- 
traction regions, vignetting and point-spread function correc- 
tions. We used the latest available spectral redistribution matrix 
(swxpc(3tol2_20010101v008.rmf). 
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The spectral fitting was performed using XSPEC0 version 
12.3, adopting an absorbed power-law model. We find the fol- 
lowing best-fit 0^ ed = 1.00 for 359 de grees of freedom) param- 
eterfl absorption Nn = (5.7 ± 0.3) x 10 21 cirT 2 and photon 
index r = 1.78 + 0.05 (see Fig.[2j». No evidence for emission or 
absorption lines was found by inspecting the residuals from the 
best-fit models (see Fig.|2j. The addition of a blackbody compo- 
nent is not required. We used the resulting averaged flux (cor- 
rected for the absorption) of ~2.2 x 10" 11 erg crrT 2 s _I to derive 



the conversion factor 1 count s 



1.1 x 10" 10 erg cm" 2 s" 1 



in 

the 0.3-10 keV energy range. The mean count rate of each ob- 
servation is reported in Tableland the values are plotted against 
the orbital radio phase in Fig. [3] First we provide an orbit-by- 
orbit view, then all the data are plotted together and finally 
they are averaged over a ~2.5 day bin corresponding to 1/10 
of the radio orbital phase. The source luminosity ranges from 
~3.4 x 10 33 of 2 erg s to ~1.6 x lO 34 ^ erg s _1 , where we indi- 
cate with the distance in units of N kpc. 




1 2 
Energy (keV) 

Fig. 2. Swift fXRT cumulative spectrum of LSI +61° 303 (see 
Sect. |2] for details). Top: data and best-fit power-law model. 
Bottom: residuals from the best-fit model in units of standard 
deviation. 



3. Discussion 

LS I +61° 303 is one of the most intriguing Be/X-ray binary sys- 
tems because of its periodic radio and X-ray emission and of 
its strong gamma-ray emission. Two main models have been 
proposed to account for the peculiarities of LSI +61° 303 and 
of the similar bina ry L S 5039 (orbital period P or b = 3.9 days; 
iParedes et ai1l2000l) ; see lMirabell d2006l) for a recent review. The 
first model assumes that the systems are powered by accretion on 
either a neutron star or a black-hole, and that non-thermal parti- 
cles are accelerated in relativistic jets. The second supposes that 
LSI +61° 303 and LS5039 are similar to the rotation-powered 
system PSRB1259-63, composed by a millisecond radio pul- 
sar orbiting around a Be star (P or b = 1237 days and pulsar spin 
period of 47.7 ms; e.g.. lKaspi et alJ[T995l) . In this scenario, the 
non-thermal emission arises from the shock resulting from the 
interaction between the relativistic pulsar wind and the stellar 



2 See http://heasarc.gsf cnasa.gov/docs/xanadu/xspec/ . 

3 Spectral errors are given at the 90% confidence level for a single 
interesting parameter. 
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Fig. 3. Swift/XRT background-subtracted count rates folded as 
a function of the radio orbital phase (in the 0.3-10 keV energy 
range). The count rates have been converted to unabsorbed fluxes 
(in units of 10~ n erg cirT 2 s _1 ) in the same energy range using 
the conversion factor given in Section [2] Top: on each panel, the 
data, averaged over a single observation, are plotted for differ- 
ent orbits (in the fifth panel all data are plotted simultaneously). 
Bottom: all data binned at 0.1 in phase are plotted simultane- 
ously. See the online edition of the article for a color version of 
this figure. 



wind originating from the companion star. The SwiftpiKT mon- 
itoring of LS I +61° 303 over five orbital periods provides more 
elements to constrain the parameter space of the source. 

We observed flux fluctuations up to a factor of «3 from 
the mean value at intra-hour scales. An example of this vari- 
ability, already ob served with various X-ray instruments (e.g., 
Sidol i et al.l 120061) . is given in Fig.Q] This kind of variability 
has also been observed in LS 5039, and it is generally inter- 
preted ^as_be2ng_du£jo_TOriatiOTis or anisotropics in the stellar 
wind (Bos ch-Ramon et aLll2005h . The folding of the XRT count 
rates at the ~26.5 days period (Fig. strongly suggests that this 
"flickering" is superimposed on an overall orbital modulation. 
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We note that the erratic behavior of LSI +61° 303 should be 
carefully considered when assessing the orbital modulation pat- 
tern using observations of the source collected at different times. 

Previous multi- wavelength campa i gns undertaken with dif- 
ferent instruments dTavlor et al.l 119961 : Harrison et al. 2000) in- 
dicate that the X-ray emission peak tends to precede the ra- 
dio outburst by a few days. However, the exact phase of the 
X-ray maximum was difficult to estimate because of the poor 
sampling of the orbit of LS I +61° 303. Recently, folding more 
than 600 INTEGRAL pointings performed discontinuously with 
IBIS/ISGRI from March 2003 to July 2006 for a total exposure 
time of Ms, and Ros s iXTE/ ASM data from a similar time- 
span, Hermse n & Kuiperl d2007l) found a maximum at phase 
~0.55 in the 1-100 keV band (in a greement with the results 
found by iParedes et ail d 19971) and iGrundstrom et all d2007l) 
based on RossiXTE/ASM data only). 

By binning the folded light curves in different ways, we es- 
timate that a broad X-ray peak occurs around phase 0.65 (see 
Fig. [3]), while the radio outburst is ex pected at phase a;0.8- 
0.9 dGreeorv et alii 19991: lGregorvll2002l) . The ma ximum X-ray 
flux is not detected at periastron (<p 0.23; ICasares et alj 
120051). where the accretion rate should reach its maximum value 
( Mart i &Paredeslll995l) . 

In the framework of the model invoking a shocked pulsar 
wind, X-ray flux peak and radio outburst could be in principle 
ascribed to the same particle population. Althou gh significant 
inverse Compton losses should be at work (e.g.. lDubusll2006t 
ISidoli et al1l2006l) . it is worth noting that the delay of «0.2 in 
phase (~5 x 10 5 s) between X-ray and radio peaks is compati- 
ble with synchrotron-dominated cooling rate of electrons with 
typical Lorentz factor of y ~ 10 5 assumi ng a magn etic field 
of ~0.1 G in the post-shock distribution. Dubus (2006) predicts 
higher values (a few Gauss) for the magnetic field at a stand-off 
distance of R s ~ 10 11 cm from the pulsar, but their theoreti- 
cal estimate relies on as yet poorly constrained parameters such 
as pulsar spin-down energy and the ratio of magnetic to kinetic 
energy in the pulsar wind. In this context, the fact that both X- 
ray and radio peaks are not detected at periastron is not surpris- 
ing, considering that shock strength could strongly depend on 
the anisotropic geometry of the pulsar wind, and not only on the 
separation between the two stars. 

The data collected with the XRT detector can be described 
by a featureless power-law spectrum with photon index T ~ 1.8. 
Similar photon indices were also measured with other X-ray 
telescopes ("e.g. ISidoli et alJ [20061) . This value fits reasonably 
well in a shock acc eleration scenari o where electrons and 
positrons cool slowly dChevalierll2000h . while spectra of high- 
mass accret ion-powered X-ra y binaries are generally harder 
(r * 1; e.g JWhite etafll 19831) . 

We remark that the spectral characteristics of LS I +61° 303 
are similar to those of the Be/X-ray binary A 0538-66 
dWhite & Carpenter! 1 19781) in its present low-luminosity state. 
This system, situated in the Large Magellanic Cloud, was dis- 
covered as a source of X-ray flares showi ng a recurrenc e of 16.65 
days, interpreted as the orbital period (Skinne 3 119801) . Optical 
obser vations also unveile d a super-orbital modulation of ~421 
days dAlcock et alll200l|). A0538-66 contains a 69 ms X-ray 
pulsar (ISkinner et al.lll982h and it showed iron line -emission in 
at least one low-amplitude flare dCorbet et alJll997h . 

In recent years A 0538-66 entered a quiescent phase, with 
outbursts 2-3 orders of magnitude below the e arly levels and low 
persistent luminosity (Campana et al. 200 2]). Using wXMM- 
Newton observation performed in 2002. lKre tschmar et al. ( 2004) 
found a spectrum well fitted by a featureless power-law model 



with r = 1.9 + 0.3 and estimated a 0.3-10 keV luminosity of 
(5 - 8) x 10 33 erg s _1 . The mechanism presently powering the 
faint X-ray emission of A 0538-66 is unclear. It could be ei- 
ther a low-rate accretion or the pulsar rotational energy through 
the shock between the pulsar wind and the inflowing mate- 
rial, as proposed for LS I +61° 303. To our knowledge, no data 
concerning the radio or GeV-TeV emission of A 0538-66 are 
available to further investigate the similarities and the possible 
connection between this source and the "gamma-ray binaries" 
LS I +61° 303, LS 5039, and PSRB1259-63. 
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